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Abstract 

The present view about the structure of measons is that they are a quark-antiquark 
system. The mass spectrum corresponding to this system should, in principle, be given by 
chromodynamics, but this turns out to be a complex affair. Thus it is of some interest to 
consider relativistic systems of particle-antiparticle, with a simple type of interaction, which 
could give some insight on the spectra we can expect for mesons. In the present paper we 
carry this analysis when the interaction is of the Dirac oscillator type. We show that the 
Dirac equation of the antiparticle can be obtained from that of the particle by just changing 
the frequency u into -u;. Following a procedure suggested by Barut we derive the equation 
for the particle-antiparticle system and solve it by a perturbation procedure. We thus obtain 
explicit expressions for the square of the mass spectra and discuss its implications in the 
meson case. 


1 Introduction and summary 

It is well known 1 * that mesons are considered as formed by a quark-antiquark systems where, in 
many cases, 1 * the particle and antiparticle are of the same mass i.e. mi, ud, du, dd; ss etc. 

The calculation of the mass spectra of mesons within the framework of quantum chromody- 
namics would be a complex affair 2 *. Thus it is of some interest to consider relativistic systems of 
particle-antiparticle, with a simple type of interaction, that could give us some insight in the type 
of spectra that we can expect for mesons. 

In the present paper we intend to carry this analysis when the interaction is of a Dirac oscillator 
type 3,4 *. We begin in section 2 by considering the positive and negative energy solutions of the 
one particle Dirac oscillator problem 3,4 *, and show that the equation for the anti-particle can be 
derived from that of the particle if we change the frequency w of the oscillator to — u>. 

In section 3 we consider the Barut 5 * procedure for deriving a single equation for n-free rel- 
ativistic particles of spin 1/2, and generalize it to n particles with Dirac oscillator interactions 
with different frequencies u/ f , s = 1,2, ...n. We then apply it to the particle-antiparticle case where 
n = 2 and a?! = — = u/. 
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In section 4 we reduce our equation, which has four conqjonents, to just a single one, and 
proceed to show how to solve the ratter by perturbation theory. ; . . , _ 

In section 5 we derive explicitly the square of the mass spectra of our particle-antiparticle 
system to first order perturbation theory, and proceed to draw the square of the mass level 
scheme as function of the total angular momentum of the parity (-1) J or -(-1) J as well as of 
the number of quanta N of the oscillator, for different values of u. 

Finally, in the concluding section, we also give the square of the mass spectra of the mesons 
and show that, while quite different from our present theoretical analysis, it could, as in the three 
quark case of baryons 6,7 *, give a better agreement if other interactions are also considered. 


2 The Dirac oscillator equation for a particle and for an 
antiparticle 

The single particle Dirac oscillator equation was suggested by the replacement 3 * 


p -» p - iwx/3, 


in the Dirac free particle expression 8 * giving j-ise to 


i(dtp/dx°) = a • (p - iux/3) -I- /?j 


where i° is the time and u the frequency of the oscillator, all in the units 

h = m = c = I 

where m is the mass of the particle and c the velocity of light. Note furthermore that 

/0 ff\ „ (l 0 \ 


-C o) ■ *-«•,). 


(2.4a, 6) 


where <r is the Pauli spin vector. v WkkzlS z ‘ : .i -t? 

We require now the solutions of (2.2) both for positive and negative energy where E will denote 
its absolute value. For positive energy we can write 

^ jivn/ — Z * 1 ( *?_ .t'I 




where i/>f are the large and small components depending only on the coordinates and, from 
(2.4a,b), we obtain 


(E - l)rpt = <* ■ (P + tux) 

( E + 1 )rl>% = a • (p - twx) ikt, 


(2.6a) 

( 2 . 66 ) 


so that from the second equation 

= (E + l) _1 [<r • (p - iwx)]^ , 
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and substituting in the first we get 


where 


( E 2 — 1)0* — [p 2 + wV — 3w — 4u;L • S]0*, 


( 2 . 8 ) 


L = x x p , S = a/2. 


(2.9a, 6) 


Clearly then 0* is given by 

0i + = \n(i, ^);m) = 0 = R N( {r) [{^rn -a,^<T|>m)J^ m _ <r (tf,^)xi/ 2 <r], 


( 2 . 10 ) 


while to obtain 0j we have to carry out the operation (2.7). In (2.10) Rm{r) is the radial 
function given in terms of Laguerre polynomials, < | > a Clebsch-Gordan coefficient, V/ m _„(0, y>) 
a spherical harmonic, and Xi/2» a spin 1/2 state with projection <j. The N indicates the total 
number of quanta while t,j = l ± j are the orbital and total angular momentum. 

Now we turn our attention to the negative energy states where we can take 




from which 


-(£+i)0r = 

(-£ + 1 ) 0 ^ = 


O • (p + io/x) 
a • (p — tu>x) 


02 > 

07 > 


so that 


while xl ?2 satisfies 


07 = -(£ + 1) _1 \a • (p + iu?x)j 07 , 
(E 2 — 1)07 = p 2 + w 2 x 2 + 3 cj + 4wL ■ sj 07, 


( 2 . 11 ) 

(2.12a) 

(2.126) 

(2.13) 

(2.14) 


so it is again given by the ket \N(t, \)jm) = 0 of (2.10). 

The particle state 0 P is the positive energy one 0* which from (2.5) can be written as 10) 


0p = 0 + = 


0 

(£ + 1) *[a- (p - twx)]0 


exp(— »£x°), 


(2.15) 


with 0 given by (2.10). 

For the antiparticle state 0„ we follow Bjorken and Drell 9 * by taking the conjugate of 0“ of 
(2.11) and apply to it 

o 1 )’ (2I6) 
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0a = 1720 ’ = 


exp(— iEx°) 


thus getting 10 ’ 


iff 20* 

( E + 1) -1 <7 • (p + twx)(icr 2 0*) 


(2.17) 


as 

( — it 72)<J p ’( — ttfj) = ° > ( — *>2)( t<r 2) = (2.18) 

Furthermore as 

*> 2 X 11 = -Xi-J > t<y 2 Xi-i=Xii (2.19a, b) 

we see that 

r 1 i 

i&2$ m = ^ft/w( r ) XI |{^ m “ a ' 2 <T ^ 771 ^”l) m+3 -a 

<r= — 1/2 L 

( 2 . 20 ) 

so changing ff,m into — <7, -m and using properties of the Clebsch-Gordan coefficients 11 ’ we obtain 

t<720'=(-l) m+i -V (2-21) 

Thus, (except for the phase factor (-l) m+;_ ' which is irrelevant) the state 0 a of the antiparticle 
is the solution of the Dirac oscillator equation (2.2) when we change u> by — uj. 

In the next section we consider a Poincare invariant equation for the two body system of 
particle and antiparticle. 


3 Equation for the particle-antiparticle system with Dirac 
oscillator interaction 

As in previous publications*' 10 ’ we start from the Dirac equation for n-free particles 


X>,P , + /3,)0 = £0, (3-1) 

j=i 

where E is the total energy for the system, a,,/3, are direct products such as 

& = 7®/.../®/?®/...®/, (3-2) 

with /3 in s tl * position, while p, is the momentum of the s th particle. 

Following the analysis of Barut 5 ’, we showed that the Poincare invariant form of the equation 

(3.1) is*- 10 ’ 

Efr.(7^ + i)lv- = o. (33) 

We first explain all the symbols appearing in (3.3). The index /i takes now the values \i — 

0,1, 2, 3, and „ . 

7 ° = = 1,2,3, (3.4a, b) 


i 

i 

1 


i 



408 



with 0, a, given by ( 2 . 4 a,b). tor n particles we have 7 = 1 , 2 , ...n, given by direct products 
such as (3.2). 

We also introduce the concept 5 ' of unit time like four vector (u M ) = (u 0 , u\, u 2 , U3), which 
means that in some reference frame it can tkae the value (1,0, 0,0). With the help of u M we define 
the Lorentz scalars 


r = nw “<•). 

r=l 

r f = (7X)- , r, 


(3.5a) 

(3.56) 


where repeated indeces are summed over fi = 0, 1 , 2, 3. Note that in (3.5b) just eliminates 

the corresponding term in T of (3.5a) so T s is still in product form. 

The terms in our equation (3.3) are then fully defined and we proceed now to look at it in the 
frame of reference where (u M ) = ( 1 , 0 , 0 , 0 ), where it takes the form 


where 


(r° £ po, + ^2 r°(7, • p, + 

^ 3=1 3=1 



r° = IK 

r® = (7, 0 )- 1 


= 0®0®....®0 = B, 

r° = 0 ® /?..../? ®i® /?.... ® 0, 


Multiplying (3.6) by T° and using 0 2 = 7 , 7 , = /?a, and (T 0 ) 2 = / we obtain 


- P° + £>. • p, + A)U = 0, 

• 3=1 


where we put the time like component P 0 of the four momentum 


= p M i + p M 2 + ■••• + Pun » 


(3.6) 


(3.7a) 

(3.76) 


(3.8) 


(3.9) 


in its contravariant form — P° as our metric tensor is 


9^ = 0 if n ± u, g n = g 22 = 933 = ~9oo = 1. 


(3.10) 


Clearly we then recover equation (3.1) if we interpret P° as E, as is usually done. 

The Barut equation (3.3) will then provide the starting point for the one involving n-particles 
with Dirac oscillator interactions of frequency w,,s = 1 , 2 , ...n. To proceed in this direction we 
could replace p M# ; p = 0 , 1 , 2 , 3; s = 1 , 2 , ...n in ( 3 . 3 ) by a linear function of p M , and x„, as was done 
in (2.1) for the one particle problem. We note though that while p M , commutes with the total four 
momentum P v of (3.9), which is a genereator of the Poincare group, x M , does not. Thus it is more 
suggestive 4 ' to use the translationally invariant coordinates x'^, defined by 

x' M , = x„.-X„, ( 3 . 11 a) 
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where 


Xp = n _1 (x Ml + x M 2 + .... + x Mn ). (3.116) 

Furthermore we would like that the resulting equation, in the frame of reference where the 
center of mass is at rest i.e. P, = 0, i = 1,2,3, should depend only on a single time as is the case 
in Eq. (3.8). Thus it is convenient to use the transverse coordinates x^ defined by 

K. = <■ - (f , ‘K.)f , Jp,p r r' n , ( 3 . 12 ) 

which has the property that in the center of mass frame, where Pi = 0, i = 1, 2, 3, x' 0a = 0. 

With the above restrictions we could obtain from (3.3) a Poincare invariant equation with a 
Dirac oscillator interaction if we make in it the replacement 


Pm -*■ Pm ~ 


K. r > 


(3.13) 


where T is defined as in (3.5a), and we assignate a different frequency u>, for each particle. We 
then arrive at the equation 


n 


E 



- iw.x^n + 1 1^ = 0, 


(3.14) 


where we make the choice 4 ^ for the unit time like vector in T and T 9 of (3*5) as 

u u = P^-PvPT 112 (3.15) 

In the center of mass frame of reference, where P, = 0, i = 1, 2,3, we have {u„ } = (1000) and 
so by a reasoning similar to the one that leads from (3.3) to (3.8) we obtain the equation 


{ - -P 0 + E [<*' ■ (Pi - + 0s] = o, 


(3.16) 


where B is given by (3.7a) while a,,j3, are direct products of the form (3.2) and 


Pi = P. -n P, 


(3.17) 


becomes identical to p, in the center of mass frame. 

As P° is the total energy of the system, the rest of the expression (3.16) is then the mass 
operator 4 *, which we will designate by M, for the n particles interacting with Dirac oscillators 
of frequencies u,,s = 1,2, ...n. If we are dealing with the particle-antiparticle system n = 2 and, 
from the discussion of the previous section u\ = — = w, so we get 4 * 

M = (l/\/2){(ori — ttj) • p — + a 2 ) - x ]#) + (0i + &) (3.18) 

where 4 * 

p = (l/\/2)(pi - p 2 ), x = (l/\/2)(xi - x 2 ), (3.19o, 6) 
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-c ;m: :> 

—a :x 7 > 

»-« ’,)•(! !) 
*-(! !)•(! i) 


0 

*1, 

0 

0 \ 

*1 

0 

0 

0 

0 

0 

0 


0 

0 

*1 

0 / 


0 0 *2 0 \ 

0 0 0 <r 2 I 

<r 2 0 0 0 I ’ 

^ 0 0 2 o 0 / 


(I 0 0 0 \ 

0-/00 
”00/0’ 
\0 0 0 -// 

// 0 0 0 \ 
0/0 0 
“00-/ 0 ’ 
\0 0 0 -// 


(3.20a) 


(3.206) 


(3.21a) 


(3.216) 


// 0 0 0 

B -(I 0 \ (I 0 \ _ 0 — / 0 0 

B ~ Vo -i) ® Vo -// o o -/ o 

\0 0 0 / 

and the wave function can be written as 4 * 


0 


4 Solution of the eq. (3.18) by a perturbative procedure 

Denoting by /i the eigenvalue of the mass operator of (3.18) and making use of (3.19-3.23) we 
obtain the equation 

0 <T\ • (p + itJX) — (T 2 * (p - iwx) 0 \ /0n\ 

o x • (p — iux) 0 0 —02 • (P + lwx ) 02i 

• (p + iux) 0 0 01 (p- *«x) 0ia 

0 — <7j • (p — iwx) ffi • (p + iux) 0 /V t />22 / 

<(/i-2) 0ii > 

0 ^21 (4.1) 

M 012 ’ 

\ (A* + 2) 022 / 





411 



where p = — xV. 

Introducing now the creation and annihilation operators by the definitions 
tj = (1/\/2)(w 1/2 x - iuT 1/2 p),( = (1/\/2)(lj 1/2 x + iu~ 1/2 p), 
we get the equations 

■ 1/2 (<T\ *1 *2 ■ { \ f 02 i \ f (t* - 2) i \ 

\*2'{ C\ • Tf ) \ 012 / V (a* + 2) 022 y’ 

-iw 1/2 f * 1 * 2 f 01,N ) = u(^ 21 ) 

\*2 »? < T \ ■ ( ) \022 / ^ V 012 / 


(4.2a, 6) 


(4.3a) 

(4.36) 


Multiplying (4.3a) by \x and substituting in it (4.3b), we obtain, after some straightforward algebra, 
that 

- 2u 0 


U) 


’ A 

D' 

’011 | _ 

D 

A. 

. 022 J 


0 


f + 2/i 


][:;;]• 


where 


while 


A = 2tf ■ ( + 3 - L • (<r t - <r 2 ), 

D = 2(S-t,) 3 + 2(S-() 2 -(i-t,)-(f.(), 

L = x x p = — 1(*7 x (),S = ^(<r, +tr 2 ), 
and extensive use was made of the relation between Pauli spin matrices i.e. 

<7,<jj = 6,j + ie tJ k( 7k- 

It is convenient to substitute 0n,022 by 0 + ,0_, through the relation 

m U-i !][':]■ 


011 
L 022 


so that equation (4.4) becomes 

[A-D 

w l 0 


0 

A + D 


ntj- u; 7 m 


(4.4) 


(4.5a) 

(4.56) 

(4.6a, 6) 


(4.7) 


(4.8) 


(4.9) 


Writing the two equations in 0 + , 0_ explicitly and eliminating 0_ between them we obtain for 
0+, which from now on we denote simply by 0, the equation 

— (4 + 2Au)fi 2 + w 2 (A 2 — D 2 — AD + Z);4)j 0 = 0 (4.10) 

Unfortunately, because of the term D of (4.5b), this problem is not exactly soluble as was 
the case of two particles, i.e. ui\ = = u>, discussed in reference 4. We note though that the 

operator in (4.10) contains the frequency u> as a parameter. As this frequency is given in units 


! 

I 


I 
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of the rest mass 1 of the particle, we expect w to be small as compared with 1 , as is the case in 
nuclear physics. We can then begin by disregarding the term in u> 2 and so our equation becomes 


M 2 - (4 + 2 Au) <j> = 0, 


(4.11) 


so our first objective will be to find the eigenstates and eigenvalues of the operator A given by 
(4.5a). 

To achieve our purpose we start by introducing the ket 


\N{£,s)jm) = Rwt(r) Y,(tm - <T,s<j\jm)Y lm .'{9,v)x*T, (4.12) 

<7 


where all the functions and symbols are defined as in the paragraph following ( 2 . 10 ), except that 
now the spin function \ lo corresponds to s = 0 or 1 and not the 1/2 appearing in ( 2 . 10 ). 

We note from (4.5) that the operators in (4.10a,b) are invariant under reflections, i.e. change 
of 17 , £ into -t), and thus the parity of the states (4.12), which is (-1)*, is a good quantum 
number. Considering then separately the states (4.12) in which t — j, and those in which l = j ± 1 , 
we find by straightforward Racah algebra 11 * that the eigenstates of A of ( 4 . 5 a), which we shall 
denote by fo, are 


4> 0 = (l/\/2) j^|7V(j,0);m) ± \N(j, l)j'm)J for parity (— 1 ) J , (4.13a) 

4>a = | N(j ± 1 , l)j'm) for parity - (— 1 ) J , (4.136) 

The corresponding eigenvalues of p 2 , which we denote by /ij, are given by 

= 4 + 2 u ;{(2 N + 3) ± \j(j + 1 )] 1/2 |, (4.14a) 

for parity ( — l)- 7 , while for parity — ( — l)-* we get 

Mo = 4 + 2w(2JV + 3), (4.146) 

and thus we have a complete solution of the problem (4.11). 


Our interest though is in the equation (4.13) which we can solve by a perturbation procedure. 
We first define 


W = m 2 , 

(4.15a) 

H 0 = 4 + 2 Au)y 

(4.156) 

H' = u\A 2 -D 2 -AD + DA), 

(4.15c) 

so equation (4.10) becomes 


{W 2 - WH 0 + H')4> = 0. 

(4.16) 

We then, as for example in Schiff book 12 *, replace H' by A H' where A 

is a parameter and write 

W = Wo + XW, +A 2 W 2 + ..., 

(4.17a) 

4> = <t>o A A^i + A 2 ^2 + 

(4.176) 
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(4.17c) 


where W 0 = nl of (4.14) and is given by (4.13). 

From (4.17a) we obtain 

W 2 = Wq + X(2W 0 W 1 ) + \ 2 (2W 0 W 2 + W 2 ) + 


SO 


that using (4.17) we see that, up to first order in A, (4.16) takes the form 


W 0 (W 0 - H 0 )<Po 


+ A 


(2 W 0 W l - W X H 0 + H')<t > o + W 0 (W 0 - H 0 )<t>i 


+ ... — 0, 


(4.18) 


where each of the square brackets must vanish l2 b For the first one this is automatic as from 
we have 

H 0 <t>o = 


( 4 . 11 ) 

( 4 . 19 ) 


From the second square bracket, when we take its scalar product* 2 ^ with we obtain 

Wx = -WQ l (4> 0 ,H'<t> 0 ), ( 4 . 20 ) 

where we made use of the hermitian character of Ho and of Eq. (4.19). 

Thus to first order in perturbation theory, when we take, as usual l2 \ A = 1, we have that 


/i 2 = ^-Mo 2 (^o>^o) + ... 


( 4 . 21 ) 


where /i 0 is given by (4.14), <t> 0 by (4.13) and H’ is (4.15c). In the next section we calculate this 
fi 2 explicitly. 


5 Square of the mass spectra of the particle- antiparticle 
system 

To determine the square of the mass /i 2 , given to first order perturbation theory by (4.21), we 
need to calculate the scalar product (<^o, H f <i> o). As <£o has a definite number of quanta N, which 
is indicated in (4.13), we need only to consider that part of H* in (4.15c) that does not change the 
number of quanta. The terms AD , DA in (4.15) change the number of quanta by ±2, as indicated 
in (4.5), so we can disregard them. The A 2 is diagonal in the basis 4>q of (4.13) and its contribution 
to the scalar product in (4.20) is 

(A), AVo) = {(27V + 3) ± \j(j + 1)] ,/2 } 2 for parity(-iy (5.1a) 

Vo, A Vo) = (27V + 3) 2 for parity - ( -l) j ( 5l6 ) 

For the D 2 operator, where D is given by (4.5b), the only terms that contribute to its expec- 
tation value with respect to 4>q i.e. that do not change the number of quanta TV, are 

4(S ■ f*) 2 (S ■ O 2 + 4(S • 0 2 (S • »?) 2 - 2(S • « ■ i) ~ 2(S • ( ) 2 (i, ■ if) 

- 2 (t, • n)(s • O 2 - 2(f • <)(s • n i) 2 + fa • n)(t •*) + (*• £)fa • »)• ( 5 2) 
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To evaluate the matrix element (4>o, we need then to determine the matrix elemtns of 

(S-i7) 2 ,(S'{)M»r*7)i({-0 with respect to states of the form \N{£,s)jm) of (4.12). This is done 
in the Appendix using results for the matrix elements of (S • if), (S • £) given in reference 13 as well 
as the operator form of the harmonic oscillator states of reference 14. 

We finally arrive at the following results for the /x 2 of (4.21): For the states of parity ( — l) 2 ,/x 2 
can be expressed as a function of using both signs i in the of (4.13a) i.e. 

H 2 ± (N,j,u) = 4 + 2w{(2JV + 3) ± \j(j + 1)] 1/2 } 

— (u, 2 /4){[2 N + 3) ± + l) l/2 ] 2 - 2 [N(N + 3) - j(j + 1) + 3]} + ... (5.3) 

For the states of parity the two cases of orbital angular momentum £ = j + 1 or i = j - 1, 

have to be written separately. We shall distinguish them from (5.3) by putting a bar above the 
fi 2 and an index + or - when £ is respectively j + 1 or j — 1. Thus we obtain 

ft = 4 + 2w(27V + 3) - (u> 2 /4) j(2./V + 3) 2 - 2 \j(j - 1) + N(N + 3) + 1]} + ... 

(5.4a) 

fil = 4 + 2ui{2N + 3) - (w 2 /4)|(27V + 3) 2 - 2 \j(j + 3) + N(N + 3) + 3]} + ... 

(5.46) 

In all of these cases we keep only terms up to w 2 so that /x 0 2 appearing as a coefficient of the scalar 
product in (4.21) is replaced by (1/4). Note furthermore that as N = 2n + £ where £ is the orbital 
angular momentum, we see that for /x 2 {N,j,u) of (5.3) the N takes the values N = ;,;+2,;'+4, ... . 
On the other hand for /x+(JV, j,u) of (5.4a) N takes the values N = j + 1 ,j + 3, j + 5, ..., while 
for p, 2 _(N,j,(jj) of (5.4b) we have N = j - 1 ,j + 1 ,j + 3 ,j + 5,.... 

In Fig. 1 we graph /x^ of (5.3) for ui — 0.1 as function of j, indicating the sign ± to which the 
level corresponds on its left hand side while on the right hand side we give the value of N. Note 
that when j = 0 the state |^V(0, 1)00) does not exist so that the normalized <f>o of (4.13a) reduces 
to |7V(0,0)00). In this case the spin is 0 as indicated int the left hand side of the level j = 0, and 
not the mixture ± in (4.13a). 

In Fig. 2 we graph ji\ of (5.4) for u> = 0.1 as function of j. The orbital angular momentum 
( = j ± 1 is indicated on the left of the levels and the total number of quanta N is given on 
the right. Note that for j = 0, \N( — 1,1)00) does not exist so that we have only \N(l, 1)00) 
corresponding to £ — 1 indicated on the left hand side of the levels with j — 0. Also the levels 
with N = £ = j — 1 are unique as indicated in the corresponding left hand side of the levels. For 
the other states | N[j ± 1, 1 )jm) both values £ = j + 1 and j - 1 are possible, and because of the 
first order corrections, i.e. the term in u> 2 of (5.4) they are separated by uj 2 (2j + 1). In Fig. 2 we 
give on the left hand side the values £ = j + 1 and £ = j — 1 of the paired levels. 

The parity denoted by a script V is given in all three figures i.e. V = ( — 1) J or V = — ( — 1) J . 

In the next section we discuss the comparison of our reuslts with squares of the mass spectra 
for non-strange mesons 1 '. 
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6 Comparison with the meson spectra 

In page 37 of reference 1 there is Meson Summary Table in which, in the first two columns, are 
given the non-strange mesons with the code name (a single letter with an index indicating the total 
angular momentum j), mass p in MeV , j , isospin 7, and parity defined there 1 ) as P = — ( — 1)*, 
as well as a charge conjugation number C = (-1)' + *, where l is the orbital angular momentunm 
and s the total spin of the quark-antiquark system. 

From the above information we can get the square of the mass p 1 in units (GeV) 2 , the j, s, I 
as well as ( — 1)*, so that, ( = j or £ = j ± 1. As mentioned before our notation for parity will be 
the script V = ( — 1)* i.e. V = ( — 1) J ' or —( — 1)-'. 

The information given in the previous paragraph is summarized in Figs. 3 to 6, where in the 
abscissa we have j = 0,1,2, 3,4,5 and p 2 in the ordinate. We note first that for comparison 
with the particle-antiparticle system, with a Dirac oscillator interaction, we have divided the 
information according to parity V = (-l) j or -(-1) J in Figs. 3,4 or Figs. 5,6 as V is an integral 
of motion of our problem. Furthermore the isospin I — 0, 1 is completely independent from the 
Poincare group and thus of the analysis previous sections. We could then consider that in our 
theoretical p\ of (5.3) or p ± of (5.4), u> is a function of 7 and another one could be added to 
these p±,p± variables. Thus it is convenient to graph separately the levels for 7 = 1 ,"P = ( — 1) J 
(Fig. 3); 7 = 0,V = (-1) J (Fig. 4); 7 = 1, V = -(-l) 2 (Fig. 5); and 7 = Q,V = -(-1) J (Fig. 6). 

In figures 3 to 6 we put on the left hand side the name of the meson and where there are 
several of the same name we distinguish them by primes i.e. Xj, ir 1 ,,*" etc. In Figs. 3,4 we put on 
the right hand side the total spin s = 0 or 1 and furthermore we differentiate the two values by 
using either a full or dashed line. In Figs. 5,6 only the names of the mesons appear on the levels 
as the spin s is always 1. 

We now wish to make a qualitative comparison between the theoretical figures 1 to 2 and the 
experimental ones for mesons in figures 3 to 6. 

We begin with parity “P = ("l) 2 where we should compare Fig. 1 with Figs. 3,4. For angular 
momentum j = 0 the comparison in Figs. 1 and 3 is quite good and the spin is s = 0. For all 
other levels we cannot compare because in Fig. 3 they have definite spin, while in Fig. 1 we have 
50% each of admixtures of s = 0 and s = 1, in or out of phase, as shown in (4.13a). This clearly 
shows the need to add other interactions in our equation (3.14) that are Poincare invariant, which 
we shall discuss below. Note also that even for j = 0 the comparison between Fig. 1 and Fig. 4 is 
not good, even if we change the scale in Fig. 1, but this can be attributed to the fact that the tjq 
mesons may have and admixture of ss where s has different mass from u, d. 

Turning now our attention to parity V = — ( — 1) J we could compare Fig. 2 with Fig. 5 or 6. 
In Fig. 5 i.e. 7 = 1 the information is sparse though for j = 1 we see what may be a pair p',p" 
which, with a change of scale, could resemble the pair for j = \,N = 2 in Fig. 2. This is more 
noticeable when we compare, for j = 2, Fig. 2 with Fig. 6, where we see a number of pairs i.e. 
,V = 3 corresponding f' v f % ; JV = 5 to /%', f[ v \ and if we had graphed it, N = 7 to N = 9 

to / 2 v ' // ,/ 2 vm . Note also that for parity V = -(-1) J the lowest level in Figs. 5, 6 appears for 
j = 1 and not j = 0, and this is also true in the theoretical calculations of Fig. 2. 

Clearly though in all cases we would have to modify our starting hypothesis in a similar way 
as we did in the baryon example 7 ). We want to do this in a Poincare invariant way and thus it is 


convenient to introduce the four vector 7 ) 


w M = 

(6.1) 

where K 9T could be either J aT , L aT ,S 9T defined respectively by 


J 9T = L aT + S' T , 

(6.2a) 

l" - £(*>: - *.>:), 

J=1 

(6.26) 

s" = oy <)E(*;-«), 

(6.2c) 


»=l 


where xj,pj,7j; r = 0, 1,2, 3; a = 1,2 are the ones in section 3 of this paper with n = 2. 

We consider now the Poincare invariants 

W 2 = W“W„ , P 2 = — P„P M , (6.3 a, 6) 

as well as the T of (3.5a), which we combine in the form 

r(W 2 /P 2 ), (6.4) 

and in the center of mass frame t.e. when P, = 0,t = 1,2, 3, it reduces to 7 ) 

J\L\S\ (6.5a, 6, c) 


depending on whether K 9T is equal to J ar , L 9T , S 9T . 

Thus, as indicated in Eq. (2.9a) of reference 7, our equation (3.16) could be modified to 

{-P° + M + aJ 2 + bL 2 + cS 2 }rl> = 0 (6.6) 

where M. is given by (3.18) and a, 6 , c are, so far, arbitrary real constants. 

If b = c = 0, as J 2 is am integral of motion of the operator M, we have that the new mass, 
which we call ft', is related with the old one by 

fi' - aj (j + 1)] =/i 2 , (6.7) 

and thus the new mass spectra is 

fi' = n + aj(j + 1), (6.8) 

where we have at our disposal the parameter a with which we can adjust the spectra corresponding 
to different j's. 

When we have 6 , c also different from zero, as L 2 ,S 2 are not integrals of motion, our only way 
to proceed is by considering the matrix of the operator 

M + aJ 2 + bL 2 + cS 2 , (6.9) 
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with respect to the states \N(l,s)jm), where; is fixed, and diagonalizing this matrix up to certain 
maximum number of quanta N max . 

Calculations of this type, done by L u is Benet, when a = 6 = 0, i.e. only with an S 2 term, 
allow us to break the wave functions of (4.13a) into states of definite spin i.e. s = 0 or 1, but 
keeping some of the ordering as in Fig. 1, so that its more comparable with the meson spectra in 
which the spin is given. 

We do not wish though to consider the more general operator (6.9) in this paper. To begin 
with, when dealing only with the operator M of (3.18), we have a single parameter, the frequency 
u, and also as a scale in our calcualtions the mass m of the quarks. If we go to the operator (6.9) 
we have the parameters w, a, 6, c, plus the m. Thus we can of course adjust the meson spectra 
better, but it may be meaningless. 

Thus we conclude by stating that the particle-antiparticle system with a Dirac oscillator in- 
teraction, may give some insight on the meson spectra, and that is all we aspire to achieve in the 
present paper. 
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7 Appendix 


To obtain the matrix elements of (<£oi D 2 <t> o) we need first those of the factors in each product 
appearing in (5.2). Those of (S • tf) 2 can be obtained from the ones of (S • if) given in (3.20) and 
(3.24) of reference 13, and thus we have 


r i >/2 

(N + 2 (j, l)jm|(S • tl) 2 \N(j, 1 )jm) = - [(AT j + 3 )(N + 2 - j ) j 
(N + 2(j + 1, l)jm|(S • *l) 2 \N(j + 1, l)jm) = — {N + ; + 4)(AT + 1 -»] 7 [>/(2j + 1) 


(/1.1a) 


(A.16) 


ll/2r l 1 /? 

(N + 2(j - 1, l);m|(S • ti) 2 \N(j + 1, l)jm) = [(JV - j + 3 )(N + 1 - »j \j(j + 1)J (2; + l)’ 1 

(Ale) 

(N + 2(j + 1, l)jm|(S -tl) 2 \N(j - 1, l)jm) = [(JV + ;' + 2)(N + j + 4)] [j(j + 1) 


ll/2 


( 2 J + 1)- 1 
(A.ld) 


r i 1/2 r 

{N + 2(j - 1 , l)jm|(S • ri) 2 \N(j - 1 , l)jm) = - [(iV j + 2){N -j + 3)J [(j + l)/(2; + 1) 


(A.le) 


From the hermitian conjugates of (5.1) we obtain those of (S • £) 2 *- c - 


(N - 2 {j, l);m|(S • *) 2 | N(j, 1 )jm) = - [(N + j + 1)(AT - j)] 


1/2 


(A. 2a) 


{N - 2 (j + 1, l);m|(S • t) 2 \ N(j + 1, 1 )jm) = - [(AT + j + 2)(N - j - l)] 1 ^ [j/(2j + 1)] 

(N - 2 (j + 1, l)jm\(S ■ t) 2 \N(j - 1, l)jm) = [(N - j + l)(N - j - 1)] + l)]^' + l)" 1 

(/ - ’ 

(N - 2 (j - 1, l);m|(S • t) 2 \N(j + 1, l)jm) = [(JV + j)(N + j + 2)] 1/J [;(> + 1)]' *(2 j + l)' 1 

r i ,/J r 

(N - 2 (j - 1, l)jm|(S • - 1, l)jm) = - [(N j)(N - j + l)j [(j + l)/(2j + 1) 


(A. 26) 


(A. 2c) 


(A. 2d) 


(A.2e) 


Finally, from the operator form of the harmonic oscillator states, given in reference 14, we have 
that i/2 

{N + 2{e, 3 )jm\ii-il\N{t,s)jm} = - [(JV + * + 3)(JV - * + 2)] , (A.3a) 


419 


and from its hermitian conjugate we get 

[ 1 1/2 
(W + /+l)(jV-/)J . (,4.36) 

With the help of these expressions we obtain straighforwardly the n 2 ± of (5.3) and p 2 ± of (5.4). 
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9 Figure Captions 

• Fie. 1. We graph (A of (5.3) for u = 0.1 as function of j, indicating the sign ± to which 
the level corresponds on its left hand side while on the right hand side we give the value of 
N. Note that when j = 0 the state \N(0, 1)00 > does not exist so that the normalized of 
(4.13a) reduces to |JV(0,0)00 >. In this case the spin is 0 as indicated in the left hand side 
of the level j = 0, and not the mixture ± in (4.13a). 

• Fig 2 We graph ji 2 ± of (5.4) for u = 0.1 as function of ;. The orbital angular momentum 
l = j± 1 is indicated on the left of the levels and the total number of quanta N is given on the 
right. Note that for j = 0, |iV(-l, 1)00 > does not exist so that we have only \N(1 1)00 > 
that corresponds to t = 1 indicated on the left hand side of the levels with j = 0. Also the 
levels with N = £ = j — 1 are unique as indicated in the corresponding left hand side of the 
levels. For the other states \N{j ± 1, 1 )jm > both values i = j + 1 and j - 1 are possible, 
and their separation is given by w 2 (2; + 1), thus both values j — 1, j + 1 appear on the left 
hand side of those levels. 

• Fig. 3. The square of the masses of the mesons are given as functions of j for isospin 7 = 1 
and parity V = {-\ ) j . The name of the meson is given on the left hand side with an index 
j and upper primes if there are several of them. Full lines correspond to spin 0 and dashed 
to spin 1 as indicated on the right hand side. 

• Fig 4 The square of the masses of the mesons are given as function of j for isospin 7 = 0 

and parity V = (-l) 2 . The name of the meson is given on the left hand side with an index 

j and upper primes if there are several of them. Full lines correspond to spin 0 and dashed 
to spin 1 as indicated on the right hand side. 

• Fig. 5. The square of the masses of the mesons are given as function of j for isospin 7 — 1 

and parity V = — (— 1) J . The name of the meson is given on the left hand side with an index 
j and upper primes if there are several of them. The spin is always 1. 

• Fig. 6. The square of the masses of the mesons are given as function of j for isospin 7 = 0 

and parity V = — (—l)*'. The name of the meson is given on the left hand side with an index 
j and upper primes if there are several of them. The spin is always 1. 
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